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ABSTRACT 

We have measured the B hadron energy distribution in decays using a sample of 
semi-leptonic B decays recorded in the SLD experiment at SLAC. The energy of each 
tagged B hadron was reconstructed using information from the lepton and a partially- 
reconstructed charm-decay vertex. We compared the scaled energy distribution with 
several models of heavy quark fragmentation. The average scaled energy of primary B 
hadrons was found to be < xe^ >= 0.716 ± 0.011(stat.) l;o;o22 (syst.). 
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1. Introduction 



The production of heavy hadrons (H) in e"'"e~ annihilation provides a laboratory for 
the study of heavy-quark (Q) jet fragmentation. This is commonly characterised in 
terms of the observable = '^En/y/s, where is the energy of a B or D hadron 
containing a b or c quark, respectively, and ^/s is the cm. energy. In contrast to 
light-quark jet fragmentation one expects |l[ the distribution of xeh^ E>{xe„), to peak 
at an X£;^-value significantly above 0. Since the hadronisation process is intrinsically 
non-perturbative D{xeh) cannot be calculated directly using perturbative Quantum 
Chromodynamics (QCD). However, the distribution of the closely- related variable xeq 
= IEqI ^ can be calculated perturbatively ^, || and related, via model-dependent 
assumptions, to the observable quantity D{xeh)\ ^ number of such models of heavy 
quark fragmentation have been proposed 0. Measurements of D{xe„) thus serve 
to constrain both perturbative QCD and the model predictions. Furthermore, the 
measurement of D{xeh) different cm. energies can be used to test QCD evolution, 
and comparison of D{xeb) with D{xeo) '^^^ ^e used to test heavy quark symmetry P). 
Finally, the uncertainty on the forms of D{xeo) ^^"^ D{xeb) iiiust be taken into account 
in studies of the production and decay of heavy quarks, see eg. more accurate 
measurements of these forms will allow increased precision in tests of the electroweak 
heavy-quark sector. 

Here we consider measurement of the B hadron scaled energy distribution D{xeb) 
in Z° decays. Earlier studies |10| used the momentum spectrum of the lepton from 
semi-leptonic B decays to constrain the mean value < xeb > and found it to be 
approximately 0.70; this is in agreement with the results of similar studies at -^i = 29 
and 35 GeV [|ri|]. In more recent analyses |J2|, |1^ the scaled energy distribution D{xeb) 
has been measured by reconstructing B hadrons via their B D/X decay mode; we 
have applied a similar technique. We used the precise SLD tracking system to select 
jets containing a B — D/X decay, where the charmed hadron D was identified semi- 
inclusively from a secondary decay vertex formed from charged tracks. Each hadronic 
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vertex was then associated with a lepton Z (/ = e or /x) with large momentum transverse 
to the jet direction. Neutral energy depositions measured in the hermetic calorimeter, 
as well as the energies of charged tracks, that were not associated with the DZ system 
were subtracted from the jet energy to yield the reconstructed B hadron energy. This 
measurement technique may be useful to B-lifetime or B-mixing analyses where 
the proper time t = Lj \/^^ — 1, where 7 = EbIvtlb^ vhb is the B hadron mass and 
L is the decay length, must be known accurately. We then compared the B energy 
distribution with the perturbative QCD and phenomenological model predictions. 



2. Apparatus and Hadronic Event Selection 

The e^e~ annihilation events produced at the resonance by the SLAC Linear Col- 
lider (SLC) were recorded using the SLC Large Detector (SLD). A general description 
of the SLD can be found elsewhere |jl5|. This analysis used charged tracks measured in 
the Central Drift Chamber (CDC) and in the Vertex Detector (VXD) [|17|, energy 
clusters measured in the Liquid Argon Calorimeter (LAC) |TB[, and muons measured 
in the Warm Iron Calorimeter (WIC) Electron identification utilizes CDC tracks 



and LAC clusters pO|] . 



Momentum measurement is provided by a uniform axial magnetic field of 0.6 T. 
The CDC and VXD give a momentum resolution of (Jp^/p^ = 0.01 © 0.0026p^, where 
p± is the track momentum transverse to the beam axis in GeV/c. Including the 
uncertainty on the primary interaction point (IP), the resolution on the charged- 
track impact parameter (d) projected in the plane perpendicular to the beamline is 
ad = ll©70/(p_Lv^sin^) /im, where 9 is the polar angle with respect to the beam- 
line. This results in a mean resolution on reconstructed 2-prong vertices (Section 
3) of cr^^Kj^) = 400 (25) /im for the projection on an axis along (perpendicular to) 
the vertex flight direction. The LAC electromagnetic energy scale was calibrated 
from the measured 7r° — > 77 signal the electromagnetic energy resolution 
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is ge/E ^ 0.15/^E(GeV). 

The trigger and initial selection of hadronic events are described in A set of cuts 
was applied to the data to select well-measured tracks and events well-contained within 
the detector acceptance. Charged tracks were required to have a distance of closest 
approach transverse to the beam axis within 5 cm, and within 10 cm along the axis from 
the measured interaction point, as well as | cos^| < 0.80, and p± > 0.15 GeV/c. Events 



were required to have a minimum of seven such tracks, a thrust axis polar angle 9t 
within I cos6't'| < 0.70, and a charged visible energy E^s of at least 20 GeV, which was 
calculated from the selected tracks assigned the charged pion mass. From our 1993-95 
data sample 108650 events passed these cuts. The efficiency for selecting hadronic 
events satisfying the | cos 9t\ cut was estimated to be above 96%. The background in 
the selected event sample was estimated to be 0.1 ± 0.1%, dominated by — >• r+r^ 
events. 

Calorimeter clusters used in the subsequent jet-finding analysis (Section 4) were 
required to comprise at least two calorimeter towers, each containing an energy of 
at least 100 MeV, and to have a total energy greater than 250 MeV. Electromagnetic 
clusters used in the non-B-associated neutral energy measurement were further required 
to have less than the smaller of, 25% of their energy and 600 MeV, in the hadronic 
section of the LAC. 

The efficiency for reconstructing B hadrons, the background in the selected sample, 
and the resolution of the method were evaluated (Sections 3 and 4) using a detailed 
Monte Carlo (MC) simulation. The JETSET 7.4 event generator was used, with 



parameter values tuned to hadronic e+e annihilation data p6|, combined with a 



simulation of B-decays tuned to T(45') data |2^ and a simulation of the SLD based 
on GEANT 3.21 Inclusive distributions of single particle and event topology 

observables in hadronic events were found to be well-described by the simulation ||29|| . 
There is now evidence that roughly 21% of all promptly-produced B hadrons in 
— > bb events are B** mesons since JETSET does not produce B** mesons we 
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have corrected the simulation to account for them. Using an event weighting technique 
we produced a generator-level distribution of B hadron energies in which the energy 
Eb of 20.7% of all B hadrons was adjusted to be Eb — E^^, where the pion energy 
Et^ was produced according to an isotropic 2-body decay distribution for B** Bvr^, 
assuming a B** mass of 5.7 GeV/c^. Uncertainties in this simulation of B** production 
were taken into account in the systematic errors (Section 7). 



3. B Hadron Selection 

Hadronic events were required to contain a lepton candidate within the barrel tracking 
system with |cos^^| < 0.7. We then applied the JADE jet-finding algorithm to 



the LAC clusters in each selected event to define a jet topology. With a jet-resolution 
criterion of yc = 0.07, 82.9% of the events were classified as 2-jet-like and 17.1% 
as 3-jet-like. Kinematic information based on this topological classification was used 
subsequently (Section 4) in the calculation of the B hadron energy. Events in which 
the lepton had a transverse momentum w.r.t. its jet axis, pt, of at least 1 GeV/c were 
retained for further analysis. In jets containing more than one such lepton only the 
highest-pt lepton was labelled for association with a D vertex and any lower-momentum 
leptons were used in the D-vertex-finding. 

In each selected jet we then searched for a secondary D vertex among the non-lepton 
tracks. Tracks were required to comprise at least 40 CDC hits and one VXD hit, to be 
well contained within the CDC with | cos^| < 0.70, to have momentum in the range 
0.15 < p < 55 GeV/c, and to have a transverse impact parameter, normalised by its 
error, of d/ad > 1. Tracks from and A° decays and 7 conversions were suppressed 
by requiring the distance of closest approach to the IP in the planes both perpendicular 
to, and containing, the beamline to be less than 1 cm. Two-prong vertices were first 
formed from all pairs of tracks whose distance-of-closest-approach was less than 0.012 
cm and whose fit to a vertex satisfied < 5. A multi-prong D-vertex candidate was 
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then defined to comprise the tracks in all accepted two- prong vertices in the jet, and to 
be located at the position of the two-prong vertex containing the track with the largest 
normalised transverse impact parameter d/ad- 

The tracks in each D vertex were each assigned the charged pion mass and were 
then combined by adding their four- vectors to obtain the vertex invariant mass, m^), 
and the vertex momentum vector. The vertex flight distance from the IP was projected 
onto the jet axis to obtain the quantity td- Events were retained if at least one jet 
contained a D vertex with 0.3 < m/) < 1.9 GeV/c^, ru > 0.05 cm, rr, normalised by its 
error larger than unity, and the distance-of-closest-approach between the lepton track 
and the extrapolated D-vertex momentum vector was less than 0.012 cm. The lepton 
and D-vertex tracks were then fitted to a common candidate B vertex. The combined 
D-vertex and lepton invariant mass, m^, and the projection of the vector between the 
B- and D-vertex positions onto the D-vcrtcx momentum vector, tb, were calculated. 
Events were selected in which ruB < 4.5 GeV/c^, > 0.025 cm, and rs normalized 
by its error was larger than unity. 

For the selected events, distributions of the number of tracks per D vertex, Nd, and 
of r£), iTiB, and are shown in Fig. 1. Also shown are the simulated distributions 
in which the contribution from selected true B — >■ D/X decays is indicated. In Fig. 2 the 
distributions of lepton transverse momentum with respect to the jet axis, pt, are shown 
for candidates passing all cuts except the requirement that pt be above 1 GeV/c; the 
simulated distributions are also shown, and the contributions from different processes 
are indicated. The final sample comprises 597 events, 293 in the muon, and 304 in the 
electron, channels. Using the simulation we estimate that the purity of this sample, 
defined to be the fraction of the tagged events whose identified leptons I are from true 
B — > DZX decays, is 69.2%; a further 18% of the selected events contain B decays with 
a cascade, punch-through or mis-identified lepton, and are still useful. The estimated 
composition of the bb events in terms of the B hadron species is shown in Table 1. 
The remaining 12.8% of the event sample comprises non-bb events. The efficiency for 
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selecting B hadron decays in the selected hadronic event sample is shown, as a function 
oi xeb^ ill Fig- 3; the overall efficiency is 1.1%. 



4. Measurement of the B Energies 

In each selected event we first defined the jet energies by using kinematic information. 
The 2-jet events were divided into two hemispheres by the plane normal to the thrust 
axis and the jet in each hemisphere was assigned the beam energy. For the 3-jct 
events we corrected the jet energies according to the angles between the jet axes, 
assuming energy and momentum conservation and massless kinematics. Labelling the 
jets arbitrarily 1, 2 and 3, and the corresponding inter-jet angles 623, 9i3 and $12, the 
corrected energy of jet 1 is given by: 

El = v^(sin^23)/(sin^i2 + sin ^23 + sin ^13), (1) 

with corresponding expressions for jets 2 and 3. This procedure results in improved jet 
energy resolution. 

We then proceeded to reconstruct the B hadron energy E^'^: 

rpRec 771 771 /r)\ 

— ^jet — -t^frag, {■^) 

where Ejet is the energy of the jet containing the candidate B vertex and Efrag is the 
energy in the same jet that is not attributed to the B, 

77' fchg Tpchg , rneu rnneu /n\ 

Efrag = / E f^^g + / E f^^g (3j 

where -E/^ag ^frag mcasurcd charged and neutral energy components re- 

spectively, and Z'^'*^ and /"'^^ are correction factors described below. We define E'^'^^g 
to be the sum of the energy, using the momentum and assuming the pion mass, of 
all the charged tracks in the jet excluding the candidate B- vertex tracks; -B/rog i^ 
fined to be the sum of the energy of the electromagnetic calorimeter clusters in the jet 
that are not associated with charged tracks. A cluster was defined as unassociated if 



it had no charged track extrapolating to it to within an angle Aaa from its centroid, 



where ad = y '^d + '^d ^'^'^ ^ci ^"^^ ^ti measured cluster widths in polar- and 

azimuthal-angle, respectively. The distributions of and E'^^^g ^'^^ shown in Fig. 4. 

This procedure will a priori misassign the energy of any unassociated neutral par- 
ticle from the D decay to the non-B energy Ef^ag- Similarly, the energy of any charged 
track from the D decay that is not associated with the reconstructed D vertex will be 
misassigned to Efrag- We have used our MC simulation to study these effects and show 
in Fig. 5 the correlation between the reconstructed and true values of -E/^ag Ej'^^g. 
As expected, both the charged and neutral components are typically slightly overes- 
timated by the reconstruction method. We fitted an ad hoc second-order polynomial 
to each correlation to determine an average energy-dependent correction factor, f'^^ 
(/"^^") (Eq. 3), which we applied to the non-B charged (neutral) energy component 
Ej^lg (Ej^^g) of each tagged jet in the data sample. Uncertainties in these corrections 
were included in the systematic errors (Section 7). 

We have used our simulation to estimate the resolution of the method for recon- 
structing the B hadron energy. We compared the reconstructed scaled B energy 
with the input scaled energy a;^^^ and show in Fig. 6 the distribution of the quantity 
(^x^rue _ x'^£^)/x%^^. The resolution may be characterised by a parametrisation com- 
prising the sum of two Gaussian distributions. The result of such a fit, in which the 
Gaussian centers, normalisations and widths were allowed to vary, is shown in Fig. 6. 
The narrower Gaussian of width a — 0.10 represents 65% of the fitted area, and the 
wider Gaussian of width a = 0.33 represents the remainder. It can be seen from Fig. 6 
that the population corresponding to the 'inner core' is somewhat underestimated by 
this technique since the parametrisation does not describe the central bin. We repeated 
this exercise in subset regions of x%'^^ and found the inner core resolution (population) 
to be 0.27 (84%) for 0.0 < x%-^^ < 0.6, 0.09 (70%) for 0.6 < x%l^ < 0.8, and 0.06 
(79%) for 0.9 < a;^'^^ < 1.0; as expected the resolution is better for more energetic B 
hadrons. Choosing the bin width to be roughly half of our mean resolution we show 
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the measured distribution of x^£^, D'^°'^°'{x'^£^), in Fig. 7. Also shown in this figure is 
the simulated distribution in which the background contribution from non-bb events 
is indicated. 



5. Comparison with Model Predictions 

It is interesting to compare our measured B hadron energy distribution with the the- 
oretical predictions. The event generator used in our simulation is based on a per- 
turbative QCD 'parton shower' for production of quarks and gluons, together with 
the phenomenological Peterson function ^ (Table 2) to account for the fragmenta- 
tion of b and c quarks into B and D hadrons, respectively, within the iterative Lund 



string hadronisation mechanism this simulation yields a generator-level primary 
B-hadron energy distribution with < > = 0.693Q. It is apparent (Fig. 7) that this 
simulation does not reproduce the data well; the fo^' the comparison is 36.7 for 15 
bins. 

We have also considered alternative forms of the fragmentation function based on 
the phenomenological model of the Lund group 0, the perturbative QCD calcula- 
tions of Braaten et al. ||, (BCFY) and of Nason et al. (NCM), as well as ad 
hoc parametrisations based on a function used by the ALEPH Collaboration |T2[ and 



on a third-order polynomial. These functions are listed in Table 2. 

In order to make a consistent comparison of each function with the data we adopted 
the following procedure. Starting values of the arbitrary parameters were assigned and 
the corresponding distribution of scaled primary B hadron energies, D*^'"(x^^^), was 
reproduced in our MC-generated bb event sample, before simulation of the detector, 
by weighting events accordingly. The resulting distribution, after simulation of the 
detector, application of the analysis cuts and background subtraction, of reconstructed 
B hadron energies, D^"^*" (x^^), was then compared with the background-subtracted 
* We used a value of the Peterson function parameter tb = 0.006 [52|. 
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data distribution and the value was calculated. This process was iterated to find 
the minimum in x^, yielding a parameter set that gives an optimal description of the 
reconstructed data by the input fragmentation function. This procedure was applied 
for each function listed in Table 2. The fitted parameters and minimum values 
are listed in Table 3, and the corresponding D^'^{x^£^) are compared with the data in 
Fig. 8. Each function reproduces the data. We conclude that, within our resolution and 
with our current data sample, we are unable to distinguish between these functions. 
It should be noted, however, that the optimal third-order polynomial function has a 
small negative minimum point in the region around x%^'^ = 0.2; since this behaviour is 
unphysical we did not consider this function further in the analysis. 



6. Correction of the B Energy Distribution 

In order to compare our results with those from other experiments it is necessary 
to correct the reconstructed scaled B hadron energy distribution D'^°'^"'{x^£^) for the 
effects of non-B backgrounds, detector acceptance, event selection and analysis bias, 
and initial-state radiation, as well as for bin-to-bin migration effects caused by the 
finite resolution of the detector and the analysis technique. We also corrected for the 
effects of B** decays (Section 2) to derive the primary B hadron energy distribution. 
We applied a 15 x 15 matrix unfolding procedure to to obtain an estimate 

of the true distribution 

D'^'%x%:^) = e-\x%:^) ■ E{x%l\ xTb) ■ {D^'^'^xTJ - SixTJ) (4) 

where S* is a vector representing the background contribution, E is a matrix to correct 
for bin-to-bin migrations, and e is a vector representing the efficiency for selecting true 
B hadron decays for the analysis. 

The matrices S, E and e were calculated from our MC simulation; the elements 
of e are shown in Fig. 3. The matrix E incorporates a convolution of the input frag- 
mentation function with the resolution of the detector. We used in turn the Peterson, 
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Lund, BCFY, NCM and ALEPH functions, with the optimised parameters hsted in 
Table 3, to produce both a generator-level input primary B hadron energy distribution 
j-)MC ^^true^ ^ and a reconstructed distribution as discussed in the previous 

section. In each case E was evaluated by examining the population migrations of true 
B hadrons between bins of the input scaled B energy, x^^*^, and the reconstructed 
scaled B energy, x^^. 

The data were then unfolded according to Eq. (4) to yield D'^*'°'{x%'^), which is 
shown for each input fragmentation function in Fig. 9. It can be seen that the shapes 
q£ jjdata ^^true^ differ Systematically among the assumed input fragmentation functions. 
These difference were used to assign systematic errors, as discussed in the next section. 



7. Systematic Errors 

We have considered sources of systematic uncertainty that potentially affect our mea- 
surement of the B-hadron energy distribution. These may be divided into uncertainties 
in modelling the detector and uncertainties on experimental measurements serving as 
input parameters to the underlying physics modelling. For these studies our standard 
simulation, employing the Peterson fragmentation function, was used. 

The uncertainty on the correction of the non-B neutral jet energy component -E/^ag 
(Section 4) was estimated by changing the LAC cluster-energy selection requirement 
from 100 to 200 MeV, and by varying the LAC electromagnetic energy scale within our 
estimated uncertainty of ±2.2% of its nominal value |21|. In each case the difference 



in results relative to our standard procedure was taken as the systematic uncertainty. 
A large source of detector modelling uncertainty was found to relate to knowledge 
of the charged tracking efficiency of the detector, which we varied by our estimated 
uncertainty of ±2.4%. In addition, in each bin of a;^^, we varied the estimated con- 
tribution from fake leptons in the data sample (Fig. 2) by ±25%. These uncertainties 
were assumed to be uncorrelated and were added in quadrature to obtain the detector 
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modelling uncertainty in each bin of xeb- 

As a cross-check we also varied the event selection requirements. The thrust-axis 
containment cut was varied in the range 0.65 < | cos^^t| < 0.70, the minimum number 
of charged tracks required was increased from 7 to 8, and the total charged-track energy 
requirement was increased from 20 to 22 GeV. In each case results consistent with the 
standard selection were obtained. As a further cross-check on jet axis modelling we 
systematically varied Uc in the range 0.01 < i/c < 0.15 and repeated the analysis; results 
consistent with the standard analysis were obtained. 

A large number of measured quantities relating to the production and decay of 
charm and bottom hadrons are used as input to our simulation. In bb events we have 
considered the uncertainties on: the branching fraction for bb; the rates of 

production of B^j, and B^ mesons, and B baryons; the rate of production of B** 
mesons, and the B** mass; the branching ratios for B ^ D* and B D**; the lifetimes 
of B mesons and baryons; and the average charged multiplicity of B hadron decays. 
In cc events we have considered the uncertainties on: the branching fraction for 
—>■ cc; the charmed hadron fragmentation function; the rates of production of D'^, 
and Dg mesons, and charmed baryons; and the charged multiplicity of charmed hadron 
decays. We have also considered the rate of production of ss in the jet fragmentation 
process, and the production of secondary bb and cc from gluon splitting. The world- 
average values 1^, ^ of these quantities used in our simulation, as well as the respective 
uncertainties, are listed in Table 4. 

The variation of each quantity within its uncertainty was produced in turn in our 
simulated event sample using an event weighting technique The matrices S and 
E (Section 6) were then reevaluated using the simulated events, and the data were 
recorrected. In each case the deviation w.r.t. the standard corrected result was taken 
as a separate systematic error. These uncertainties were conservatively assumed to 
be uncorrelated and were added in quadrature to obtain a total physics modelling 
uncertainty in each bin of xeb- 



12 



The model-dependence of the unfolding procedure was estimated by considering the 
envelope of the unfolded results illustrated in Fig. 9. In each bin of we calculated 
the average value of the five unfolded results, as well as the r.m.s. deviation. The 
average value was taken as our central value in each bin, and the r.m.s. value was 
assigned as the respective unfolding uncertainty. 

8. Summary and Conclusions 

We have used the precise SLD tracking system to reconstruct the energies of B hadrons 
in e+e" Z° events via the B DiX decay mode. We estimate our resolution on the 
B energy to be about 10% for roughly 65% of the reconstructed decays. The distribu- 
tion of reconstructed scaled B hadron energy, L>(x^^), was compared with perturbative 
QCD and phenomenological model predictions; the calculations of Braatcn, Cheung 
and Yuan and of Nason, Colangelo and Mele are consistent with our data, as are the 
phenomenological models of Peterson et al. and of the Lund group. The distribution 
was then corrected for bin-to-bin migrations caused by the resolution of the method 
and for selection efficiency, as well as for the effects of B** production, to derive the 
energy distribution of primary B hadrons produced by decays. Systematic uncer- 
tainties in the correction were considered. The final corrected XEg distribution D[xeb) 
is listed in Table 5 and shown in Fig. 10; the statistical, experimental systematic, and 
unfolding uncertainties are indicated separately. 

It is conventional to evaluate the mean of this distribution, < xeb >■ For ©ach of 
the five functions used to correct the data we evaluated < > from the distribution 
that corresponds to the optimised parameters; these are listed in Table 3. We took the 
average of the five values of < > as our central result, and defined the unfolding 
uncertainty to be the r.m.s. deviation. We hst in Table 4 the errors on < xeb > 
resulting from the study of detector and physics modeUing described in Section 7. We 
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obtained: 



0.716 ± 0.011(stat.) to^ll (exp. syst.) ± 0.019 (unfolding) 



where the systematic error is the sum in quadrature of the individual contributions 
listed in Table 4. It can be seen that < XEg > is relatively insensitive to the variety of 
allowed forms of the shape of the fragmentation function D{xEg)- 

Our results are in agreement with a previous measurement of the shape of the 
primary B hadron energy distribution at the resonance |]T2[, as well as with mea- 
surements of the shape |jl3| and mean value [0 of the distribution for weakly-decaying 



B hadrons, after taking account of our estimate that the latter < xe^ > value is about 
0.015 lower. Combining all systematic errors in quadrature we obtain < xe^ > = 
0.716 ± 0.011 (stat.) (syst.). 
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B species 


C (%) 


e(%) 


B„ 


43 


92 


B, 


43 


87 


B. 


10 


89 


B baryons 


4 


87 



Table 1: The composition C of true B D/X decays in the final sample; e is the 
fraction of each species whose D vertices are correctly reconstructed. In all cases the 
MC statistical errors are less than 2%. 



Function Name 


Functional form D{x) 


Reference 


Peterson 


i(l-i-^)- 

X 1—X ' 






§ 


Lund 


i(l - x)''exp(- 


'hm\/x) 




m 


BCFY 


{\-{\-r)xf 


/i(r) + x2/2(r) 


-x^h{r) + x^fi{r)] 


i 


NCM 


Jdyg{x,y)y'^{l 








ALEPH 


l+b{l-x)/^ c 


d \-2 




Hi 




x ^ X 


l-x) 






3»''^_order Polynomial 


1 + bx + cx"^ + dx^ 







Table 2: Fragmentation functions used in comparison with the data. For the BCFY 
function /i(r) = 3(3 - 4r), /2(r) = 12 - 23r + 26r2, f^{r-) = (1 - r)(9 - llr + 12r2), 
and /4(r) = 3(1 - rf{l - r + r^). 
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Function 


xVd.o.f. 


parameters 




Peterson 


14.0/11 


ej, = 0.034 ±0.006* 


0.717 


Lund 


9.6/10 


a = 1.7 ±0.2 
6 = 0.19 ±0.01 


0.743 


BCFY 


22.4/11 


r = 0.20 ±0.02 


0.705 


NCM 


15.9/11 


a = 9±2 
/? = 44 ± 8 


0.687 


ALEPH 


9.7/9 


6 = 0.0 ±1.0 
c = 0.78 ± 0.05 
d = 0.042 ± 0.004 


0.730 


3'"'^-order polynomial 


14.9/9 


b = -7.53 ±0.04 
c= 16.49 ±0.07 
d = -9.98 ± 0.07 





Table 3: Results of optimisation of fragmentation functions to the reconstructed B 
hadron energy distribution. For the NCM fit the QCD parameters were fixed at A/ 
= 200 MeV and fi = mb = 4.5 GeV. * This value of Sb refers to the B-hadron energy 
distribution; it should not be confused with the value of used as input in the JETSET 
model at the b-quark fragmentation level (Section 5), which is significantly lower. 
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Error source 


Variation 


Error (%) 


DETECTOR MODELLING 






Neutral fragmentation energy: 






cluster energy scale 


±2.2% 


+0.12 
-0.27 


min. clus energy 


lOO^Joo MeV 


+0.00 
-0.21 


Tracking inefficiency 


2.4 T 2.4% 


+0.2 
-1.0 


Lepton mis-ID background 


±25% 


+0.66 
-0.65 


PHYSICS MODELLING 






B meson / baryon lifetime 


1.55 ±0.05 / 1.10 ±0.08 ps 


+0.11 
-0.12 


B** production 


20.7 ±7% 


+0.68 
-0.10 


B** mass 


5.704 ±0.020 GeV 


+0.03 
-0.00 


/* = r(B^D*)/r(B^D) 


J -f'/3 


+0.32 
-0.00 


/** = r(B^D**)/r(B^D) 




+0.32 
-0.21 


Bn, Bd / Bs / b-baryon production 


40.1 ± 20.0% / 11.6 ± 8.0% / 7.0 ± 4.0% 


+0.51 
-0.48 


B„, Bd, Bs, b-baryon decay modes 


±lo- 


+0.11 
-0.12 


B-decay charged multiplicity 


5.3±0.2 tracks 


+0.25 
-0.16 


c- fragmentation: < xeo ^ 


0.484 ±0.008 


±0.01 


D° / D+ / Ds / c-baryon production 


56.0 ± 5.3% / 23.0 ± 3.7% / 12.0 ± 7.0% / 8.9 ± 0.5% 


±0.01 


D decay multiplicity 


Ref. H 


+0.04 
-0.05 


ss production 


±10% 


+0.37 
-0.40 


Rb 


0.2216 ±0.0010 


+0.00 
-0.01 


Rc 


0.16 ±0.01 


+0.02 
-0.04 


g — !■ bb splitting 


±50% 


+0.23 
-0.30 


g — > cc splitting 


±50% 


+0.22 

-0.25 


Total 




+ 1.32 

-1.48 



Table 4: Systematic errors on (xeb)- 
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1 / rrdrr / c\ t o 

A- 1 KJ \JLkJ I \XtAj 




kJ you. \j±. X \j±. 


TTnfolrliTiP' iTnpprf./^i'nf.v 


037 














nun 

yj . ± ±\j 


104 


041 


055 


041 


n 183 


n 105 


050 


068 


035 


n 256 


158 


076 


095 


043 


0.329 


0.248 


0.099 


0.102 


0.064 


0.402 


0.358 


0.115 


0.096 


0.074 


0.475 


0.560 


0.136 


0.095 


0.061 


0.548 


0.951 


0.167 


0.126 


0.033 


0.621 


1.489 


0.204 


0.137 


0.088 


0.694 


2.136 


0.242 


0.164 


0.171 


0.767 


3.011 


0.278 


0.164 


0.191 


0.840 


2.944 


0.285 


0.251 


0.112 


0.913 


1.460 


0.211 


0.319 


0.144 


0.986 


0.164 


0.067 


0.118 


0.041 



Table 5: The fully-corrected scaled B hadron energy distribution. 
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Figure Captions 

Figure 1: Candidate D-vertex distributions: (a) number of tracks per vertex; (b) 
vertex mass; (c) projection of the vertex flight distance from the IP along the jet axis. 
Candidate B-vertex distributions: (d) vertex mass; (e) projection along the D-vertex 
momentum vector of the vector between the D vertex and the B vertex. Data (points 
with error bars) and simulation (solid histogram); the dashed histogram shows the 
simulated contribution from true B — > DZX decays. In (a) all cuts were applied. In 
(b)-(e) all cuts were applied except those on the quantity shown, and these latter cut 
positions (sec text) are indicated by arrows. 

Figure 2: Distribution of (a) electron and (b) muon transverse momentum w.r.t. the 
jet axis in jets containing a selected D vertex and respective lepton. Data (points with 
error bars) and simulation (histogram) . The composition of the simulated distributions 
in terms of leptons from B — > I decays, cascade B — > C — > Z decays, wrongly- assigned 
leptons, promptly produced C — > 1 decays, and fake leptons is indicated. 

Figure 3: The efficiency e for selecting B hadron decays, as a function of scaled energy 
xeb- Note that the first bin (no point shown) is beneath the kinematic limit for xeb- 

Figure 4: Distribution of non-B-associated (a) charged and (b) neutral energy in jets 
containing a candidate B D/X decay. Data (points with error bars) and simulation 
(histogram). 

Figure 5: Simulated correlation between the true and reconstructed values of the non- 
B-associated (a) neutral and (b) charged energy in jets containing a candidate B — > DZX 
decay. In each bin of reconstructed energy the error bar represents the corresponding 
r.m.s. deviation in the true energy. Each line represents a fit to the correlation (see 
text) . 



25 



Figure 6: Distribution of the normalised difference between tfie true and reconstructed 
B hadron energies in simulated events. The sohd hne is a fit of the sum of two Gaussian 
distributions (see text). The two component Gaussian distributions are indicated by 
the dashed lines. 

Figure 7: The distribution of reconstructed scaled energies for B hadron candidates; 
data (points with error bars) and simulation (solid histogram). Also shown (dashed 
histogram) is the simulated contribution from non-bb events. 

Figure 8: The background-subtracted distribution of reconstructed scaled B hadron 
energy. The data (points with error bars) are compared with simulations based on six 
different input B fragmentation functions (see text) represented by lines joining entries 
at the bin centers. 

Figure 9: Data distribution of scaled B hadron energy corrected using simulations based 
on different input B fragmentation functions (see text): (a) ALEPH, (b) Peterson, (c) 
Lund, (d) BCFY and (e) NCM functions. Statistical error bars are shown; these are 
highly correlated between bins and among the five sets of results, (f) The five optimised 
functional forms used in the correction. 

Figure 10: The final corrected distribution of scaled B hadron energies. In each bin the 
statistical error is indicated by the innermost error bar, the quadrature sum of statis- 
tical and experimental systematic errors by the middle error bar, and the quadrature 
sum of statistical, experimental systematic and unfolding errors by the outermost error 
bar. Note that the first bin (no point shown) is beneath the kinematic limit for xeb- 
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